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ABSTRACT 

Context. Some high mass X-ray binaries (HMXB) have been recently confirmed as 7-ray sources by ground based Cherenkov 
telescopes. In this work, we discuss the 7-ray emission from X-ray transient sources formed by a Be star and a highly magnetized 
neutron star. This kind of systems can produce variable hadronic 7-ray emission through the mechanism proposed by Cheng 
and Ruderman, where a proton beam accelerated in the pulsar magnetosphere impacts the transient accretion disk. We choose 
as case of study the best known system of this class: A0535+26. 

Aims. We aim at making quantitative predictions about the very high-energy radiation generated in Be-X ray binary systems 
with strongly magnetized neutron stars. 

Methods. We study the gamma-ray emission generated during a major X-ray outburst of a HMXB adopting for the model the 
parameters of A0535+26. The emerging photon signal from the disk is determined by the grammage of the disk that modulates 
the optical depth. The electromagnetic cascades initiated by photons absorbed in the disk are explored, making use of the 
so-called "Approximation A" to solve the cascade equations. Very high energy photons induce Inverse Compton cascades in 
the photon field of the massive star. We implemented Monte Carlo simulations of these cascades, in order to estimate the 
characteristics of the resulting spectrum. 

Results. TeV emission should be detectable by Cherenkov telescopes during a major X-ray outburst of a binary formed by 
a Be star and a highly magnetized neutron star. The 7-ray light curve is found to evolve in anti-correlation with the X-ray signal. 

Key words. X-ray: binaries-gamma-rays: theory-stars: individual: A0535+26 



1. Introduction 

Galactic transient and variable X-ray sources have long 
been proposed as 7-ray emitters (see e.g. Rudak & 
Meszaros 1991, Cheng et al. 1991, Anchordoqui et al. 
2003). In fact, some sources observed by the Energetic 
Gamma Ray Experiment Telescope (EGRET) have al- 
ready been suggested to be associated with high mass X- 
ray binaries (HMXBs): A0535+26 (Romero et al. 2001), 
LS I +61°303 (Kniffen et al. 1997), Cyg X-3 (Mori et 
al. 1997), LS 5039 (Paredes et al. 2000). Moreover, very 
high energy 7-ray emission was recently detected from 
the microquasars LS 5039 (Aharonian et al. 2005a) and 
LS I +61°303 (Albert et al. 2006), and from the colliding 
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wind Be-pulsar system PSR B1259-63 (Aharonian et al. 
2005b). 

In accreting HMXBs the matter is usually transferred 
onto the compact object via the stellar wind of the mas- 
sive star. In addition to its regular mass loss, the primary, 
commonly of Be spectral type, has recurrent phases of vi- 
olent ejection. This ejected material arrives to the vicinity 
of the compact object giving rise to enhanced X-ray emis- 
sion (an X-ray outburst). In some cases, the HMXB con- 
tains a strongly magnetized neutron star (NS), being the 
magnetic field intense enough as to disrupt the accretion 
disk. The accreting material is then channeled through the 
held lines onto the surface of the NS (e.g. Gosh & Lamb, 
1979). The Cheng- Ruderman mechanism acts in the mag- 
netosphere of such an accreting NS to accelerate protons 
up to relativistic (multi-TeV) energies (see a description 
of the general picture in Orellana & Romero 2005, and 
references therein). These energetic protons move along 
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the magnetic field lines to impact into the accretion disk. 
Their interaction with the material of the disk initiates 
hadronic and electromagnetic showers. 

In order to make quantitative predictions about the 
high-energy emission originated in this particular kind 
of systems we have adopted as case of study A0535+26, 
which is the best studied among the transient Be/X-ray 
binaries. A major outburst of this system was observed 
in 1994 with the instrument BATSE of the Compton 
Gamma- Ray Observatory (Finger et al. 1996). Its light 
curve allows to estimate the history of the neutron star 
mass accretion rate, which is used in this work to calcu- 
late the evolution of the parameters that determine the 
high energy emission. The power injected inside the disk 
by 7-rays resulting from neutral pion decays can be then 
estimated. We show that the cascades triggered in the disk 
are dominated by Bremsstrahlung cooling of the electrons. 
In a broad energy range the development of these cascades 
can be estimated by solving the cascade equations apply- 
ing an analytical approach, the so-called Approximation A 
(Rossi & Greissen, 1941). The gamma-ray signal emerging 
from the other side of the disk is formed by the unabsorbed 
photons and by the cascade recycled radiation. 

The propagation of the gamma-rays in the next stage 
is determined by their interaction with the radiation field 
of the donor star. Detailed analysis of the anisotropic 
cascades inside massive binary systems have been per- 
formed by Bednarek (1997, 2000, 2006) and Sierpowska & 
Bednarek (2005). We have simulated the Inverse Compton 
(IC) cascades under some simplifying assumptions. For 
this purpose we have implemented a code based on the 
scheme presented by Protheroe (1986) and Protheroe et 
al. (1992), that allows us to estimate the emerging spec- 
trum at high energies. Detectability of the TeV flux by 
the new generation of Cherenkov telescopes can be then 
discussed. 

The structure of the paper is as follows: first, we ex- 
amine how 7-rays can be produced in transient Be/X-ray 
binaries; Section 2 presents a brief review of the Cheng- 
Ruderman magnetosphcric model and the gamma-ray 
generation in such a context. Section 3 provides some in- 
formation on the archetypal system A0535+26. In Section 
4 we present the results of our time-dependent calculations 
adopting the input parameters determined during a ma- 
jor outburst in this source. Sections 5 and 6 deal with the 
cascades induced in the disk by the relativistic protons, 
and by the emerging 7-rays in the ambient photosphere. 
We close in Section 7 with a discussion. 

2. Gamma-ray emission from HMXB containing a 
highly magnetized neutron star 

2.1. The Cheng-Ruderman magnetosphere model 

Cheng & Ruderman (1989, 1991) have studied the mag- 
netosphere of a neutron star surrounded by a Keplerian 
accretion disk. When the inner parts of the disk ro- 
tate faster than the star, inertial effects lead to a mag- 



netospheric charge separation. A gap entirely empty of 
plasma is formed and it separates oppositely charged re- 
gions. In this gap E-B ^ and a strong potential drop 
AV is established. For strong shielding conditions in the 
gap the potential drop can reach the static gap value 

AV~4x 10 14 /T 5/2 ^ 1,7 R-^L% 7 B^ 7 V, where 
M* ~ 1.4 M is the mass of the neutron star, Rq is its 
radius in units of 10 6 cm, B12 is the magnetic field in units 
of 10 12 G, and L37 is the X-ray luminosity in units of 10 37 
erg s _1 . The parameter (3 = 2R^/Rp L ~ 1 is twice the ra- 
tio of the inner accretion disk radius to the Alfven radius. 
At a radius r ~ Rq the disk is disrupted by the magnetic 
pressure and the matter is channeled by the magnetic field 
lines to the neutron star polar cap, where it impacts pro- 
ducing hard X-ray emission. Cheng & Ruderman (1989, 
1991) have analyzed some possible configurations of the 
static gap for the aligned rotator, but a realistic three- 
dimensional model for the misaligned case (ft x ji ^ 0) 
has not been developed yet. For illustrative purposes, a 
sketch of the mechanism is presented in Figure [TJ 

For typical parameters in Be/X-ray binaries, AV^ ~ 
10 13-14 V (Romero et al., 2001). Protons entering into the 
gap from the stellar co-rotating region will be accelerated 
to multi-TeV energies and directed toward the accretion 
disk along the magnetic field lines. The maximum current 
that can flow through the gap can be determined from the 
requirement that the azimuthal magnetic field induced by 
the current does not exceed that of the initial field B: 
>/max ~ B±RIRq 2 , where £?* is the magnetic field and i?* 
is the radius of the neutron star (Cheng & Ruderman, 
1989). 

Note that the impact of energetic protons on the disk 
does not occur only at a singular point but over a cer- 
tain area that is radially extended, following the spread 
of the field lines that pervade the gap and are dragged 
by the disk. Some authors incorporate a beaming factor 
(AS1/47T « 0.3) in the escaping photon luminosity to take 
into account in phenomenological way of the anisotropy of 
the emerging 7-ray emission (Cheng et al. 1991, Romero 
et al. 1991). 

The accretion disk is partially threaded by the mag- 
netic field lines of the NS. The screening effect of the cur- 
rents induced in the disk surface can be taken into account 
through a screening factor n < 1. Then, the magnetic field 
is B z = rjn/r 3 inside the disk, where fi is the magnetic 
dipolar moment of the NS. An expression for the location 
of the inner radius Rq where the disk is disrupted by the 
magnetic forces, as a function of the screening factor, was 
found by Wang (1996). 

An azimuthal component of the magnetic field is orig- 
inated by the relative rotation, and it is expected that 
B z ~ £?0 close to Rq. In Orellana & Romero (2005) two 
situations were considered concerning the screening of the 
disk i.e. 77 = 1 and 0.2. The way the results of calculations 
are affected by this change was discussed in that paper. 
Throughout the present work we consider r) = 0.1, for a 
partially screened disk. 
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Fig. 1. Cheng- Ruderman mechanism at work in the mag- 
netosphere of a HMXB containing a strongly magnetized 
neutron star (adapted from Romero et al. 2001) 

2.2. production 

Relativistic protons entering into the accretion disk in- 
teract with target protons to produce secondary charged 
and neutral 7r-mesons. The subsequent neutral pion decays 
lead to gamma-ray production. Charged pions of multi- 
TeV energies are able to interact with protons of the disk 
before the decay time (Anchordoqui et al. 2003). Such in- 
teraction occurs with a cross section a^± p ~ 22 mb, in 
this energy range. 

The differential 7-ray emissivity is 



q 7 (E 7 ) = 2 



<2V° {En) 



(1) 



where S™ n (S 7 ) = E 1 + ^ 



Applying the (5-function ap- 



proximation for the differential cross sectiorQ(Aharonian 
& Atoyan, 2000), the pion emissivity becomes 

q 7T o{E 7I ) = 4-7T 8(Ek - nE kin )J p (E p ) a pp {E p ) dE p 

= ^J p (m p c 2 + ^)a pp (m p c 2 + ^) (2) 

for proton energies greater than the energy threshold 
Eth = 1.22 GeV. Here, k is the mean fraction of the 
kinetic energy -Ekin = E p — m p c 2 of the proton trans- 
ferred to a secondary meson per collision. For a broad 
energy region (GeV to TeV) k ~ 0.17. The total cross sec- 
tion of the inelastic pp collisions is well approximated by 
a ~ 30 (0.95 + 0.061n(£ kin /GeV)) mb. 

In order to compare the relative importance of the pro- 
ton energy losses in the accretion disk we have evaluated 
the cooling rates for pp interactions and synchrotron ra- 
diation. The latter process can take place because of the 
magnetic field lines are strongly twisted in the accretion 
disk. Hence, for a proton, even if initially moving along 
the field line, its inertia can force it to depart from the 
line inside the disk, where it starts to emit synchrotron 



1 This approximation considers only the most energetic neu- 
tral pion that is produced in the pp reaction aside of a fireball 
composed by a certain number of less energetic 7r-mesons of 
each flavor. See a discussion in Pfrommer and En/31in (2004). 



radiation. The cooling rates are given by (e.g. Begelman, 
Rudak & Sikora, 1990) 



^synchrj- 1 _ CTTm e B 2 E p 



6nm p c 3 



[t p p p ] = Ka pp cn p , 



(3) 



(4) 



where K ~ 0.5 is the inelasticity (fraction of the proton 
energy lost per interaction), and <tt is the Thomson cross 
section. For the particular case studied in the next sections 
we have that t pp <C i^ynchr over a wide range of energy 
(see Figure [5]) . In addition, it should be noted that the 
crossing time t CIOSS ~ 2/i/c of a proton traversing the disk 
is much longer than t pp , so the protons effectively cool by 
pp interactions in the disk, being the synchrotron losses 
negligible. 

Since the protons arrive to the disk with basically the 
same energy acquired at the gap, estimated as E p = e AV, 
the proton injected spectrum can be then estimated as 



T (p\ **™X(F F\ protonS 



(5) 



Performing the integral in expression ([I]) we obtain 
2 T 

^ "max 



enR 2 , 



&pp(E p ^) 



ph 



^k 2 (E p - m p c 2 ) 2 - m^c 4 atom cm 3 eV s ' 



(6) 



for E 1 between E x = 0.5^(1 - iv/c) and E 2 = 0.5^(1 + 
Utt/c) and q 1 {E 1 ) — outside this range. Here is the 
velocity of the pion with energy E n = n(E p — m p c 2 ). The 
luminosity released in the disk through 7r°-decays results: 



t-2h 

L°(E 7 ) = ttRI I q 7 (E^)E^ 2 n p 
Jo 



(7) 



where n p is the number density of the protons in the 
disk, and X pp is the mean free path of a proton enter- 
ing into the disk. Taking into account the opacity of the 
disk (absorption in the Coulomb field of nuclei, for which 
<7 7P ~ 10~ 26 cm) we can estimate the luminosity that es- 
capes. Energetic photons that are absorbed initiate elec- 
tromagnetic cascades inside the disk, which are discussed 
in Section [6l The unabsorbed high-energy photons con- 
serve their original spectrum, completely dominated by 
the most energetic photons, as can be noted from ([6]) and 



max 

7 



kE„ 



([7]). The maximum energy is given by E. 

The photons escape along the original arrival direction 
of the primary protons. As it was mentioned in Section 
2.1, the protons follow the field lines along the gap, which 
form different angles with the disk surface due to the dif- 
ferential dragging. This will produce some beaming of the 
emission, but not very strong, in such a way that the duty 
cycle will render a detectable source for favorable viewing 
conditions. 
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Fig. 2. Cooling rates of high-energy protons impacting the 
accretion disk. In the calculations we have adopted B = 
2.5 x 10 6 G and n p — 1.2 x 10 21 cm~ 3 for the magnetic 
field and proton density inside the disk, respectively. 

3. A case of study: A0535+26 

Among the Be/X-ray binaries, A0535+26 has been noted 
for its remarkable multi- frequency behavior. This moti- 
vated extensive observational campaigns that have made 
of it the best known system in its class, and a good sub- 
ject for the model discussed here. A complete review of 
A0535+26 is given by Giovannelli k Graziati (1992). 

A0535+26 is an eccentric system formed by a slow X- 
ray pulsar (104 s) and the BOIII-type star HDE 245770. 
The orbital period is 111 days. The system has been fre- 
quently observed to undergo transient X-ray flares with a 
wide range of intensities. In fact, very recently it under- 
went another major outburst that was detected with in- 
struments of the Rossi X-ray Timing Explorer (Coe et al., 
2006). The simultaneity of X-ray activity with transitions 
between infrared states strongly suggests that these out- 
bursts occur in relation with increments of the accretion 
rate onto the neutron star. The variations in the accretion 
rate are thought to be related to the presence of a tran- 
sient circumstellar disk around the Be donor star (Haigh, 
Coe & Fabregat, 2004). 

The distance to A0535+26 is estimated in 2.6 ± 0.4 
kpc. The spectral feature detected at X-ray energies near 
110 keV can be interpreted of cyclotron origin, indicating 
that the neutron star has an intense polar magnetic field 
of - 10 13 Gauss (Finger et al. 1996). 

Romero et al. (2001) suggested that the Cheng- 
Ruderman mechanism acting in A0535+26 could be a 
plausible explanation for the gamma-ray source 3EG 
J0542+2610. In the next sections we will develop a more 
quantitative model for the possible 7-ray emission of 
A0535+26. 

4. Time dependent model 

A major outburst of A0535+26 was detected at hard X- 
rays (20-100 keV, reaching a peak flux of ^8 crab) during 



February-March 1994, with no simultaneous observations 
at lower energies. A sequence of normal outbursts was de- 
tected before and after this giant one (see Stella, White 
& Rosner 1986, for a description of different outburst 
types). The peak of the hard X-ray flux of the 1994 major 
outburst occurred after the periastron passage. Analysis 
of the QPO features and of the inferred spin-up history 
strongly supports the idea of the formation of a transient 
accretion disk around the neutron star (Finger et al. 1996). 
The hard X-ray emission detected is thought to be orig- 
inated through the collision of the accreted matter onto 
the surface of the neutron star, at the polar caps. The 
accretion disk would emit mostly at soft X-ray energies, 
likely with a higher flux. 

We have fitted a smooth curve to the flux values re- 
ported by Finger et al. (1996). At the peak of the outburst 
(t = in our notation) the disk luminosity was assumed 
to reach L — 10 38 erg s _1 . This normalization results from 
an extrapolation of the flux detected by BATSE to lower 
energies (~ few keV). The light curve was used to esti- 
mate the evolution of the relevant physical parameters of 
the system, assuming e ~ 10% as the radiative efficiency 
The radiation pressure is assumed to dominate over other 
pressure components (magnetic field, gas). This fact de- 
termines the structure of the inner region of a standard 
accretion disk, which is assumed to be partially threaded 
by the magnetic field of the NS. 

In the panels of Figure [3] and H] we show the tempo- 
ral evolution of the inner radius of the accretion disk, the 
half-height of the disk, the z-component of the magnetic 
field, the particle density and the photon density inside 
the disk assuming black body emission, the average energy 
of such photons, the injected energy into the disk in the 
form of relativistic particles, the maximum energy of the 
first-generation gamma-rays inside the disk, the opacity 
of the disk to proton propagation, the gamma-ray lumi- 
nosity produced inside the disk as a consequence of the 
7r°-decays, the opacity of the disk to these high-energy 
gamma-rays, and the luminosity of the gamma-rays that 
finally emerge from the other side of the disk. All values 
are given at Rq. The expressions used to perform these 
calculations can be found in Orellana & Romero (2005). 
Figure [2] shows the comparison of the cooling rates given 
by and Q, adopting the highest values for B and n p 
reached during the outburst. 

The total energy injected in the disk by the proton 
bombardment is huge, as it is the power released through 
neutral-pion decays (i.e. L° ~ 10 37 erg s^ 1 ). The opti- 
cal depth of the disk to the propagation of these photons 
follows a time evolution similar to the rate of accretion of 
mass, reaching t~ (P ~ 27 at the peak of the X-ray outburst. 
Therefore at high energies a deep valley in the light curve 
is found, in anti-correlation with the X-ray behavior. 

5. Cascades inside the disk 

The protons accelerated in the gap reach the disk following 
the magnetic field lines. The particles involved in the fur- 
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Fig. 3. Temporal evolution of the physical parameters of the accretion disk, t — indicates the maximum of the X-ray 
outburst. l)Inner radius i?o, where the disk is disrupted by the magnetic forces. The following values are estimated at 
i?o- 2) Half height of the accretion disk. 3) Vertical component of the magnetic field inside the disk for the assumed 
screening factor. 4) Numerical density of protons inside the disk. 5) Numerical density of the thermal photons emitted 
by the accreting matter. 6) Characteristic energy of the mentioned photons. 



ther cascades move nearly in the same direction of the ini- 
tial protons. We consider a longitudinal treatment of the 
cascade. It starts at the point of the 7— photon production 
(after the ir° decay) and develops along the proton arrival 
direction. This direction is not well determined since we 
have estimations for B z and can argue for ~ B z in- 
side the disk, but the radial component of the magnetic 
field B r , which determines the o h s angle, remains as a 
free parameter (see Figure [1]). The latter component of 
the magnetic field inside the disk depends on the angle of 
misalignment between Cl and p, (the vector of the axis of 
rotation and the magnetic dipolar moment), and on the 
inward drag of the field lines by the conducting accreting 
matter. We follow the cascade inside the disk considering 
a pitch angle of 45 degrees with the z direction (due to 
B^ ~ B z ). 

The energy losses of the electrons that take part in the 
cascades inside the disk, result from the interaction with 
the surrounding ambient (see in Figure [3] the evolution of 
the physical parameters of the disk). By comparing the 
cooling times that characterize the rate of the relevant 
processes (Bremsstrahlung and Inverse Compton) we de- 
termine which dominates. 

The ratio between these cooling times (tsr/iic) is 
shown in Figure [5] as a function of the electron energy. 



Bremsstrahlung dominates during the entire outburst. 
Therefore we can apply some simplifying hypothesis to 
solve the cascade equations. The analytical expressions 
under Approximation A obtained by Rossi and Greisen 
(1941) for the differential spectra are implemented here. 
At energies E 3> m e c 2 Z _1 / 3 /a collision processes and 
Compton effect can be neglected and the asymptotic for- 
mulae used to describe radiation processes and pair pro- 
duction. The average behavior of the cascade under these 
assumptions can be completely solved, and the results ob- 
tained are in good agreement with more recent numer- 
ical solutions obtained by Aharonian & Plyasheshnikov 
(2003). 

The natural unit of thickness in shower theory is the 
radiation length Xq, given by 



1 



2.318^(^+1) 



ln(183 Z-^S) 
l + 0.12(Z/82) : 



nib 



(8) 



where Z is the charge of the target nucleus and n is the 
number density. The number of the injected photons of en- 
ergy Eq = iJ™ ax that will develop cascades inside the disk 
can be estimated by N in j — q 1 n p irRQ\ pp {\ — e~ 2h / X pp)E , 
resulting in ~ 9.7. x 10 34 photons of energy E ~ 830 TeV 
per second at t = 0, and AT in j ~ 3.5 x 10 34 injected pho- 
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Fig. 4. Parameter evolution associated with the effects of proton impact onto the disk, and the luminosity generated 
after the pion decay. 1) Energy power injected in the disk by the proton bombardment, 2) Maximum energy of 
the gamma-ray photons resulting from the 7r° — > 27 decay, 3) Optical depth for a proton entering into the disk, 
4) Luminosity generated through neutral pion decay, 5) Optical depth for the most energetic photons, 6) Emerging 
luminosity, estimated as L°e~ T '< p . 



tons at the beginning of the outburst (t — —11 days) with 
E Q ~ 215 TeV. 

As it can be already noted from the optical depth T^p 
in Figure [4j the cascade traverses a thicker disk at the 
maximum of the outburst than at its beginning. This fact 
is also emphasized in Figure [5] where the end of the disk 
is indicated. The cascade traverses 100 and 22.8 radiation 
lengths, respectively, for t = and t = —11 days, where 
the unit of thickness is X ~ 3.4 x 10 4 cm and 17.8 x 10 4 
cm, respectively. The difference in grammage leads to the 
differences in the resulting spectra for the outgoing 7-ray 
photons, shown in the same figure. 

On the other hand, at the beginning of the outburst 
the emitted luminosity after the cascades develop in the 
disk is E 2 dn/dE - 2 x 10 28 (E/GeV)~ - 6 erg s" 1 . This 
contribution is added to the unaffected outgoing luminos- 
ity L° sc ~ 2 x 10 23 (£;/GeV) 2 erg s -1 , which dominates at 
energies above ~ 40 GeV. 

6. Photospheric opacity: Inverse Compton 
cascades 

The opacity to propagation of 7-ray photons after their 
generation in the disk is determined by a rather com- 
plex combination of terms. At some GeVs the absorption 
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Fig. 5. Ratio of cooling times for the main interac- 
tions of electrons in the accretion disk. The dominance 
of Bremsstrahhmg allows the introduction of the simplifi- 
cations of the so-called Approximation A for the calcula- 
tion of the electromagnetic cascades inside the disk. Two 
representative instances have been considered. 

will be dominated by the X-ray photospheric components. 
They are the emission by the polar caps of the NS (hard 
X-rays) and the accretion disk (soft X-rays) . See Orellana 
& Romero (2005) for rough estimates of the optical depth 

T~/X- 

Higher energy photons will interact with the pho- 
tons emitted by the massive stellar companion (Bednarek, 
1997, 2000, 2006). The stellar wind material can have a 
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Fig. 6. Results for the averaged behavior of the cascades in the disk initiated by photons. Left: Number of photons of 
energy greater than E in a shower initiated by one photon of energy E . The thickness of the accretion disk varies with 
time, and here this is indicated for two representative extreme times of the X-ray outburst: its beginning (t = —11), 
and the moment of maximum intensity (t = 0). Right: The luminosity that emerges at the end of the disk after the 
development of the cascades, calculated for the most energetic photons injected by ir° decays (Eq = 215 TeV for 
t — — 11 and Eq = 828 TeV for t = 0). Units of time are always in days measured from the peak of the X-ray ligth 
curve, luminosities are in erg s _1 



temperature ~ 10 4 K close to the star, adding a photon 
field component that effectively absorbs 7-rays with TeV 
energies (Bosch-Ramon et al. 2006b, Orellana & Romero 
2006), though here we do not consider such a contribution. 
The pairs created by absorption can boost lower-energy 
photons to very high energy via Inverse Compton scatter- 
ing, and the absorption of these new 7-rays leads to pair 
creation, so electromagnetic cascades are developed. 

We have computed the high-energy 7-ray spectra 
formed in cascades traversing the anisotropic stellar ra- 
diation field. Monte Carlo simulations are performed after 
developing a computational code based on the scheme out- 
lined by Protheroe (1986) and Protheroe et al. (1992). The 
applied techniques enable the use of exact cross-sections 
for all the relevant processes. The IC spectra are calculated 
as in Jones (1968), modified in similar way to Bednarek 
(1997). 

We follow the cascade initiated by the injected flux 
h(E) oc E photons per second that escapes from the disk 
at the beginning of the outburst, having energies greater 
than 40 GeV (i.e. the component reprocessed by matter 
cascades is not considered here). The absorber radiation 
field has a black-body spectrum. It comes from the B-type 
star located at a distance a ~ 10 13 cm from the neutron 
star, and characterized by T e ff ~ 30000 K and i?* = 10 12 
cm. 

We introduce the effects of the finite size of the star 
and the spatial variation of the field density by consid- 
ering the geometric configuration as in Dubus (2006). At 
distances greater than ~ 7R* the star is considered as 
a point source, whose photon density profile decays with 
the square distance. For simplicity, effects of the presence 
of a magnetic field are disregarded in these simulations, 
though it may influence the pure IC cascade by changing 
the lepton direction of propagation (Bednarek, 1997). 

The distance from the source to where we follow 
the cascades is of ~ 10 times the orbital separation. 
Specifically, photons that go beyond Ax ~ 10 a are stored 
to form the outgoing spectra. Further away, the cascade 



efficiency is reduced by the dilution of the background ra- 
diation field. The emerging secondary 7-rays depend on 
the angle of their escape from the accretion disk, # b s (see 
the sketch in Fig. Q]) . This angle is measured from the z 
direction toward the B-type star in the plane that contains 
both the primary and the NS. Some resulting spectra are 
shown in Figure [3 for different values of o bs- These spec- 
tra are obtained by sorting photons into bins of width 
A(log£') = 0.15, and imposing energy conservation to ob- 
tain the position of the curve in the energy-luminosity 
plane. Note that the more energetic part of the spectra 
is not affected by the cascades, so the photons of energy 
-E" lax still dominate the total luminosity. 

According to our results, the spectra of escaping cas- 
cade 7-rays present only a slight dependence on the angle 
of propagation, at energies above ~ 10 GeV. The general 
tendency of the results is in accordance with the expecta- 
tion of a deeper absorption closer to the star. It is clear 
that the cascade effect have led to results very different 
from those expected by just the incorporation of a factor 
exp(— t). The effect of the injection of a monoenergetic 
beam of protons in the disk is the final production of a 
power-law gamma-ray emission of photon index ~ —2.2, 
between ~ 10 GeV and ~ 10 TeV. The luminosities ob- 
tained are about 6 x 10 32 erg s" 1 at 100 GeV. The spec- 
tral features shown in the Figure [7] should be detectable 
by a Cherenkov telescope like MAGIC (A0535+26 is not 
in the field of view of HESS) if combined efforts with the 
X-ray observers are successful, in the sense that a the high- 
energy observations should be made at the beginning or 
ending of a major X-ray outburst. 

7. Discussion and summary 

We revisited the scenario of the high energy gamma-ray 
emission in an accreting magnetized neutron star that 
forms part of a HMXB. The evolution of the physical pa- 
rameters during a major X-ray outburst was calculated 
for A0535+26. The results obtained point to an anti- 
correlation between the X-ray and the TeV fluxes. 
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Fig. 7. Resulting luminosity from photospheric IC cas- 
cades, for different values of the observing angle. The slope 
of the injected spectrum is indicated by the line. 



The energy power injected in the disk by the proton 
bombardment (£inj) should perturb its structure. Under 
perturbations, the standard disk model can develop ther- 
mal instabilities (Shakura & Sunyaev, 1976). If the per- 
turbation affects the radiation pressure dominated zone, 
it will swell up exponentially, with a time scale in the 
range of seconds to minutes for an accreting NS. The prob- 
lem of the stability of a standard disk is currently treated 
with numerical simulations. Recent improvements on the 
topic were made by Teresi et al. (2004) , and by Merloni & 
Nayakshin (2006). According to them, taking into account 
the stabilizing effect of radial advection, a limit cycle type 
of behavior should be expected, probably giving rise to 
QPO features in the X-ray signal. The perturbing effect 
of proton bombardment on the disk should result in an 
oscillation of the disk height with moderate amplitude, 
as those studied by Teresi et al. (2004). If this does not 
happen, we speculate that the alternative effect of pro- 
ton bombardment could be to inflate the disk enough to 
stop the mechanism responsible of the perturbation (i.e. 
by filling the accelerating gap). Without the proton injec- 
tion, the disk deflates and eventually a kind of periodic 
behavior could also be established. 

We have performed Monte Carlo simulations of the 
electromagnetic cascades initiated in the photosphere of 
the system by the high energy 7-ray photons that emerge 
unabsorbed from the disk. The outgoing spectrum is char- 
acterized by a power law extending from ~ 10 GeV up to 
~ 10 TeV. At higher energies the original spectrum is not 
affected by cascades. The calculations presented here can 
be improved through more specific geometric considera- 
tions, and the inclusion of the effect of a random magnetic 
field component. 

Luminosity levels predicted at the TeV domain of 
~ 10 33 erg s _1 , coming from galactic sources similar to 
A0535+26 at the instances of the beginning or the end- 
ing of a major X-ray outburst, should be detectable by 
the new generation of Cherenkov telescopes. The spectral 



features imprinted by the photospheric opacity could be 
detected as well. Future observation of high-energy 7-ray 
photons emitted by HMXBs with accreting magnetized 
pulsars should be planned in quick response to alerts of 
renewed X-ray activity in the sources. 
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